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Abstract 

We study the use of the multilevel Monte Carlo technique El 
in the context of the calculation of Greeks. The pathwise sensitivity analysis [5] differenti- 
ates the path evolution and reduces the payoff's smoothness. This leads to new challenges: 
the inapplicability of pathwise sensitivities to non-Lipschitz payoffs often makes the use of 
naive algorithms impossible. 

These challenges can be addressed in three different ways: payoff smoothing using condi- 
tional expectations of the payoff before maturity Q ; approximating the previous technique 
with path splitting for the final timestep [1J; using of a hybrid combination of pathwise sen- 
sitivity and the Likelihood Ratio Method [4]. We investigate the strengths and weaknesses 
of these alternatives in different multilevel Monte Carlo settings. 

1 Introduction 

In mathematical finance, Monte Carlo methods are used to compute the price of an op- 
tion by estimating the expected value E(P). P is the payoff function that depends on an 
underlying asset's scalar price S(t) which satisfies an evolution SDE of the form 

dS(t)=a(S,t)dt + b(S,t)dW t , < t < T, 5(0) given, (1) 

This is just one use of Monte Carlo in finance. In practice the prices are often quoted and 
used to calibrate our market models; the option's sensitivities to market parameters, the so- 
called Greeks, reflect the exposure to different sources of risk. Computing these is essential 
to hedge portfolios and is therefore even more important than pricing the option itself. This 
is why our research focuses on getting fast and accurate estimates of Greeks through Monte 
Carlo simulations. 
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1.1 Multilevel Monte Carlo 



Let us consider a standard Monte Carlo method using a discretisation with first order weak 
convergence (e.g. the Milstein scheme). Achieving a root-mean square error of 0(e) re- 
quires a variance of order 0(e 2 ), hence 0(s~ 2 ) independent paths. It also requires a dis- 
cretisation bias of order 0(s~ l ), thus 0(s~ l ) timesteps, giving a total computational cost 
<9(£- 3 ). 

Giles' multilevel Monte Carlo technique reduces this cost to 0(s~ 2 ) under certain con- 
ditions. The idea is to write the expected payoff with a fine discretisation using 2~ L uniform 
timesteps as a telescopic sum. Let P be the simulated payoff with a discretisation using 2 l 
uniform timesteps, 

L 

E(P L )=E(P ) + ^(Pi-Pi-i) (2) 

1=1 

We then use Monte Carlo estimators using Nj independent samples 

1 N, 

E(P / -^ 1 )-^ = ^I(^' ) -^ (0 1 ) O) 

The small corrective term P ; — P { _ j comes from the difference between a fine and a coarse 
discretisation of the same leading Brownian motion. Its magnitude depends on the strong 
convergence properties of the scheme used. Let V/ be the variance of a single sample 
Pj — P,_j . The next theorem shows that what determines the efficiency of the multilevel 
approach is the convergence speed of V/ as / — > oo. 

To ensure a better efficiency we may modify Q and use different estimators of P on 
the fine and coarse levels of 7/ as long as the telescoping sum property is respected, that is 

L 

E(P L ) = E(P ) + £ E(P/ - Pf_ x ) (4) 
i=\ 

with 

E(p/)=E(Pf). (5) 

Theorem 1.1. Let P be a function of a solution to §\§for a given Brownian path W(t); let 
Pi be the corresponding approximation using the discretisation at level I, i.e. with 2 l steps 
of width hi = 2~' T. 

If there exist independent estimators 7/ of computational complexity Cj based on Ni 
samples and there are positive constants a > \ , j8 , c\ , cj , C3 such that 



1. E(!J) 



E(Pq) ;// = 

E{P,-Pi^) ifl>0 



2. \E(Pi-P)\<cihf 

3. V(f/) <c 2 hfNf l 

4. C,<c 3 N t h; 1 

Then there is a constant C4 such that for any e <e~ l , there are values for L and N[ re- 

L 

suiting in a multilevel estimatorY = \\ Y{ with a mean-square-error MSE = E((F — E(P)) 2 ) < 

1=0 

£ 2 with a complexity C bounded by 
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( c 4 e- 2 i/p > 1 

C<{ c 4 £- 2 (log£) 2 jfp = l (6) 
( c 4 s- 2 -^ a ifO<P<l 

Proof. Seed. □ 

We usually know a thanks to the literature on weak convergence. Results in [6] give 
a = 1 for the Milstein scheme, even in the case of discontinuous payoffs. j8 is related 
to strong convergence and is practically what determines the efficiency of the multilevel 
approach. Its value depends on the payoff shape and is usually not known a priori. 



1.2 Monte Carlo Greeks 

Let us briefly recall two classic methods used to compute Greeks in a Monte Carlo setting: 
the pathwise sensitivities and the Likelihood Ratio Method. More details can be found in 
0. 

Pathwise sensitivities 

Let S = (§n)n€[0fl\ be the simulated values of the asset at the discretisation times and 
W = ffin) n e[i,N] be the corresponding set of Brownian increments. Assuming that the 
payoff P(S) is Lipschitz, we can use the chain rule and write 

dV d r f dP(S) dS(6,W) * 

N N 
where dW = [\ &Wu and p[W) = Yl Py^k) is a multivariate Normal probability density 

k=\ k=\ 

function. 

dS 

We obtain — by differentiating the discretisation of (TTl) with respect to 6 and iterating 
od 

the resulting formula.The limitation of this technique is that it requires the payoff to be 
Lipschitz and piecewise differentiable. 



Likelihood Ratio Method 

The Likelihood Ratio Method consists in writing 

V = E [P(S)] = J P(S)p(d,S)dS 



(V) 



The dependence on 6 comes through the probability density function p(6,S); assuming 
some conditions discussed in (S], we can write 



dV 
dd 



P(S) d -^BdS 



de 



, a. dlos p(S) a „ 

P(S) |^ ) P (S)dS = E 



p(§) dlogpi§) 



de 



(8) 



N 

with d5 = nd5i and p(S) = Hp(S k \S k ^) 

k=l k=\ 



The main limitation of the method is that the estimator's variance is 0(N), becoming 
infinite as we refine the discretisation. 
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1.3 Multilevel Monte Carlo Greeks 

By combining the elements of sections [TTT] and 1.2 together, we write 



dV _ dE(P) _ dE(P L ) _ dE(P ) ^^E^-^ 

1=1 



de de de de 

As in Q, we define the multilevel estimators 

M d p (i) 



de 



y o=VL^ and Y i= N r L 



£ ( dpl {i) 



(0 



de 



where ^ g 1 » are computed with the techniques presented in section 



1.2 



(9) 



(10) 



2 European call 

We first consider a Lipschitz payoff. That of the European call is 

P={S T -K) + = max(0, S T - K) 

We illustrate the techniques by computing delta (5) and vega (v), the sensitivities to the 
asset's initial value So and to its volatility o. 



2.1 Path wise sensitivities 

We consider the Black-Scholes model: the asset's evolution is modelled by a geometric 
Brownian motion dS(t) = rS(t)dt + aS(t)dW t . We use the Milstein scheme for its good 
strong convergence properties. For timesteps of width h, 



S n+l =S n -[ l+rh + cAW n + — {AW2-h) ) :=S n -D n 



(11) 



Since the payoff is Lipschitz, we can use pathwise sensitivities. The differentiation of 
equation (jTTJ gives 



dS _ dS n+ \ _ dS„ 

dSp <9So dSp " 

dSo dS n ^i dS n ~ / .... /ati/2 7 in 

0, = — -D n +S n (AW n + a(AW; - h)) 



(12) 



, do 



do 



do 



To compute Y/ we use a fine and a coarse discretisation with Nf = 2 1 and N c = 2 l ~ 1 uniform 
timesteps respectively. We denote these by the superscripts and ( i_1 ). We take Ni samples 
to compute 



i N, 
Yi = —Y 

Nik 



dP dS Nf 
dSN f dd 



(i) 



dP dS, 



N c 



(0\ 



dS Nc de 



We use the same leading Brownian motion for the fine and coarse discretisations: we 
first generate the fine Brownian increments W = (AW\, AW2, . ■ ■ ,AWjf,) and then use W c = 
(AWi + AW2, ■ ■ ■ ,AWjv>_i +AWff f ) as the coarse level's increments. 
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Estimated complexity and analysis 



Unless otherwise stated, the simulations used to illustrate this paper use the parameters 
S = 100, K = 100, r = 0.05, a = 0.20, T = 1. 

In figure [ljwe plot (V($/)) as a function of (I /hi) in a log-log plot. We then measure 
the slopes for the different estimators: this gives a numerical estimate of the parameter /3 
in theorem Combining this with the theorem, we get an estimated complexity of the 
multilevel algorithm (table [TJ. This gives the following results : 

Figure 1: Pathwise sensitivities, European call : V(Y/)(7) 



Value var with level Delta var with level 




345G7S 345673 



Vega var with level 




5 3 4 5 G 7 8 



Table 1: Pathwise sensitivities, European call : estimated complexity 



Estimator 


P 


MLMC Complexity 


Value 


^2.0 


0(£~ 2 ) 


Delta 


^0.8 


0(£~ 2 - 2 ) 


Vega 


« 1.0 


0(£- 2 \0g£ 2 ) 



Giles has shown in that j3 = 2 for the value's estimator. For Greeks, the con- 
vergence is degraded by the discontinuity of |^ = Is>k- a fraction 0(h[) of the paths 
has a final value S which is 0(h{) from the discontinuity K. For these paths, there is 
a 0(1) probability that S$ and ^ are on different sides of the strike K, implying 

Jj^§^) ( ' ^ = 0(1). Thus V(fj) = 0(hi), and J8 = 1 for the Greeks. 



dp d ^ N f 
ds Nf de 
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2.2 Pathwise sensitivities and Conditional Expectations 



We have seen that the payoff's lack of smoothness prevents the variance of Greeks' es- 
timators % from decaying quickly and limits the potential benefits of the multilevel ap- 
proach. To improve the convergence speed, we can use conditional expectations (51 . In- 
stead of simulating the whole path, we stop at the penultimate step and then for every 
fixed set W = (AWk) ke ^ N _^, we consider the full distribution of (§n\W^). With a n = 
a (S N - 1 (W) , (N - 1 )h) and b n = b (§ N - [(W),(N— l)/z), we can write 



S N (W,W N ) =S N - l (W) + a n (W)h + b n {W)W N 
We hence get a normal distribution for (§n\ W) . 



(13) 



p(S N \W) 



exp 



with 



j% = $y-i + a (Sn-i , (N - l)h) h 
(J W = b(S N . 1 ,(N-l)h)Vh 

We can thus compute E [P (5^) \W] . Using the chain rule, we get 



(14) 



M 



V = E [P(S N )} = E w [E AWn [P(S N )\W]] « - £ E P(S[ 



m=l 



(15) 



Here with the normal probability density function, <J> the normal cumulative distribu- 
tion functions, a = (1 + rh)§N-i (W) and /3 = a VhSN-i(W), we get 



E(P(S N ) | W) = P ( +{a - K) < S > fa - K 



(16) 



This expected payoff is infinitely differentiable with respect to the input parameters. We 
can apply the pathwise sensitivities technique to this smooth function at time (N —l)h. The 
multilevel estimator for the Greek is then 



1 

Yi = -Y 



do ) 



de j 



(17) 



At the fine level we use ( fT6| ) with h = hf and Wf = (AW\, AW2, . . . , A% r i) to get 
E{P(S Nf )\Wf) . We then use 



dPf 
dd 



(i) 



dStff-i dE{P(S Nf )\W f ) dE(P(S Nf )\W f ) 



dd 



+ 



dd 



(18) 



At the coarse level, directly using E(P(Sn c )\W c ) leads to an unsatisfactorily low conver- 
gence rate of V(F/). As explained in Q we use a modified estimator. The idea is to include 
the final fine Brownian increment in the computation of the expectation over the last coarse 
timestep. This guarantees that the two paths will be close to one another and helps achieve 
better variance convergence rates. 
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S still follows a simple Brownian motion with constant drift and volatility on all coarse 
steps. With W c = (AW\ + AW2, . ■ ■ , AWa/>_3 + AW^-2) an d given that the Brownian incre- 
ment on the first half of the final step is AWfy—i , we get 



p{S Nc \%,W Nf -i) 



ex P I .-2 



2a, 



(19) 



w. 



with 



A% c = %-i (flc) + a , (N c - \)h c ) h c + b (# c - l)hc) AWivy-i 



•{$ Ne -l,(N c -l)h e )y/Q2 



From this distribution we derive IE [P(S#J|W c ,AW/v._i], which leads to the same payoff 
formula as before with a c = (1 + rh c + ctAWa^— 1)5^-1 (He) and 
j8 c = G\/Ji^Sn c -i(W c ). Using it as the coarse level's payoff does not introduce any bias. 
Using the tower property we check that it satisfies condition ([5]), 



E 



, [E [P(S Nc )\W c ,AW Nf -l] \W C ] = E [P(S Nc )\W c 



Estimated complexity and analysis 

Our numerical experiments show the benefits of the conditional expectation technique on 
the European call: 

Figure 2: Pathwise sensitivities and conditional expectations, 
European call : ¥(?/)(/) 



Value var with level 



Delta var with level 





Vega var with level 
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Table 2: Pathwise sensitivities and conditional expectations, 
European call : estimated complexity 



Estimator 


P 


MLMC Complexity 


Value 


^2.0 


0(£- 2 ) 


Delta 


« 1.5 


0(£- 2 ) 


Vega 


^2.0 


0(£-' 2 ) 



A fraction 0(y/hj) of the paths arrive in the area around the strike where the conditional 
expectation — — - — - - is neither close to nor 1. In this area, its slope is 0(h, 1/ ' 2 ). The 

coarse and fine paths differ by 0(h)), we thus have 0{y/hi) difference between the coarse 
and fine Greeks' estimates. Reasoning as in O we get Y^{K^ N {...\W)) = Oih'J 1 ) for the 
Greeks' estimators. This is the convergence rate observed for 8; the higher convergence 
rate of v is not explained yet by this rough analysis and will be investigated in our future 
research. 

The main limitation of this approach is that in many situations it leads to complicated 
integral computations. Path splitting, to be discussed next, may represent a useful numerical 
approximation to this technique. 



2.3 Split pathwise sensitivities 

This technique is based on the previous one. The idea is to avoid the tricky computation of 
E [P(5;v>)|W/] and E [P(Sn c )\W c , AWfy—i] . We get numerical estimates of these values by 
"splitting" every path simulation on the final timestep. 

At the fine level: for every simulated path Wf = (AWi,AW2, . . . ,AWjv>_i), we simulate 

a set of d final increments (AW^),-^ ^ which we average to get 

E [P(S Nf )\W f ] « \ t P ^f^ AW N f )) (2°) 
" 1=1 

At the coarse level we use W c = (AWi + AW2, . . . , AWfy,_3 + AW^-2)- As before (still 
assuming a constant drift and volatility on the final coarse step), we improve the conver- 
gence rate of V(F/) by reusing AWjvy_i in our estimation of E [P(Sty.)|iy c ] . We can do 

so by constructing the final coarse increments as (AW^), e [! ^ = (AW^-i + (Awi')); e [i^ 
and using these to estimate 

1 d n 
E(P(S Nc )\W c ) = E [P(S Nc )\W c ,AW Nf -i] ~ £P(%,(W C ,AWW)) 

" (=l 

To get the Greeks, we simply compute the corresponding pathwise sensitivities. 



Estimated complexity and choice of the number of splittings 

As expected this method yields higher values of /3 than simple pathwise sensitivities: the 
convergence rates increase and tend to the rates offered by conditional expectations as d 
increases and the approximation gets more precise. 
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Figure 3: Pathwise sensitivities and path splitting, 
European call : V(f/)(/) 




Table 3: Pathwise sensitivities and path splitting, 
European call : estimated complexity 



Estimator 


d 


P 


MLMC Complexity 


Value 


10 


^2.0 


0(£- 2 ) 




500 


^2.0 


0(£- 2 ) 


Delta 


10 


« 1.0 


0(£- 2 (\0g£) 2 ) 




500 


« 1.5 


0(£- 2 ) 


Vega 


10 


« 1.6 


0(£- 2 ) 




500 


^2.0 


0(£- 2 ) 



Taking a constant number of splittings d for all levels is actually not optimal; For 
Greeks, we can write the variance of the estimator as 



1 , / / dPf 
V(F/)=-V^(E(' 



96 



+ 



1 

Wd 



%,(V( 



(/) 

dP f 
~dd 



dP c 

~de 



\w f )) 



a) 



dPA 
de J 



\w f )) 



(21) 



As explained in section 



2.2 



we have Y Wf (E(...\W f )) = 0(h 3 / 2 ) for the Greeks. We also 
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have Kfy (Y(...\Wf)) = 0(hi) for similar reasons. We optimise the variance at a fixed com- 
putational cost by choosing d such that the two terms of the sum are of similar order. Taking 

— 1/2 

d = 0(h l ) is therefore optimal. 



2.4 Vibrato Monte Carlo 



Since the previous method uses pathwise sensitivity analysis, it is not applicable when 
payoffs are discontinuous. To address this limitation, we use the Vibrato Monte Carlo 
method introduced by Giles (4[. This hybrid method combines pathwise sensitivities and 
the Likelihood Ratio Method. 



We consider again equation ( 15 I. We now use the Likelihood Ratio Method on the last 



timestep and with the notations of section [2T2] we get 

■ E 



dV 
dd 



Eawv 



,~ d(logp(S N \W)) p 

p{s n ) Te \w 



We can write p(§n\W)) as /?(/%, <%). This leads to the estimator 



dV 



1 



dd ~ Ni £■ 



{ — — ^AW N 



111= 1 



dd 



do. 



de 



-E 



AW N 



_ ' dfx, 

d{\ogp) ^ r[nn 
do. 



w 



P(S N ) 



We compute 

A(m,i) 



de 



and 



do. 



de 

r(0 



with pathwise sensitivities. 



(22) 



(23) 



With Sp l > = S N {W {m \ Awjj>), we substitute the following estimators into ((23) 



E 



AW N 



E 



AW N 



P{Sn) 
P{Sn) 



i=i 



A(m,i) 



5 a, 



#■(«) 



u 



1 



N 



( ( &^ 1/ 

1 I MiyW 

+ - 5 



V 



In a multilevel setting: at the fine level we can use ( [23] ) directly. At the coarse level, for 
the same reasons as in section 2.3 we reuse the fine brownian increments to get efficient 
estimators. We take 



W c = (AWi +AW 2 , . . . ,AW Nf -3+AW Nf - 2 ) 



\ (AH$) ieM = (AW Nf -, + (AW#)) teM 
We use the chain rule to verify that condition ([5]) is verified on the last coarse step. With 



(24) 
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the notations of equation ( fT9] ) we derive the following estimators 



E 



P /A \ d(\ogp c ) i { „ 



E 



P{S Nc ) 



do, 



= E 
1 

d 

■E 



p/A ^ d(logj7 c ) ? ( m ) Aw 



\W c {m) 



w. 



1=1 



E 



IK 



(»>) 



P(S Nc ) d -prt\wW,AW Nf 
da w c 



-l 



(=i 



A(m,i) 



( 



\ 



\w c {m) 



+ 



A(m,i) 



(m) 



(25) 



Estimated complexity 

Our numerical experiments show the following convergence rates for d = 10: 
Figure 4: Vibrato Monte Carlo, European call : Y(Yi) (I) 



Value var with level 



Delta var with level 





Vega var with level 




11 



Table 4: Vibrato Monte Carlo, European call : estimated complexity 



Estimator 


P 


MLMC Complexity 


Value 


^2.0 


0(£- 2 ) 


Delta 


« 1.5 


0(£-' Z ) 


Vega 


^2.0 


0(£-' 2 ) 



As in section 2.3 this is an approximation of the conditional expectation technique, and 
so getting the same convergence rates was expected. 



3 European digital call 



The European digital call's payoff is P = Is t >k- The discontinuity of the payoff makes the 
computation of Greeks more challenging. We cannot apply pathwise sensitivities, and so 
we use conditional expectations or Vibrato Monte Carlo. 



3.1 Pathwise sensitivities and conditional expectations 



With the same notation as in section 2.2 we compute the conditional expectations of the 
digital call's payoff. 



E(P(S Nf )\W) = ® 



a-K 



E(P(S Nc )\W c ,AW Nf - l )=^ 



a c -K 



The simulations give figure [5] and table [5] 

Table 5: Pathwise sensitivities and conditional expectations, 
digital call : estimated complexity 



Estimator 


J8 


MLMC Complexity 


Value 


w 1.4 


0(£-' 2 ) 


Delta 


^0.5 


0{£- 2 *) 


Vega 


^0.6 


0{£-' 2A ) 
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Figure 5: Pathwise sensitivities and conditional expectations, 
digital call: ¥(?,)(/) 



Value var with level Delta var with level 




3 4 5 6 7 8 3 4 5 6 7 8 



Vega var with level 




'3 4 5 6 7 8 



3.2 Vibrato Monte Carlo 

The Vibrato technique can be applied in the same way as with the European call. We get 
figure [6] and table [6j 

Table 6: Vibrato Monte Carlo, digital call : estimated complexity 



Estimator 


P 


MLMC Complexity 


Value 


« 1.3 


0(£" 2 ) 


Delta 


^0.3 


0(£- 2 - 7 ) 


Vega 


^0.5 


0(£- 2 - 5 ) 



3.3 Analysis 



The analysis presented in section 2.2 explains why we expected j8 = 3/2 for the value's 
estimator. 

A fraction 0(y/h) of all paths arrive in the area around the payoff where 
(dK(P(§N)\W)/d§N-i) is not close to ; there its derivative is 0(hj l ) and we have 
\§N f — Sn c \ = 0(hi). For these paths, we thus have 0(1) difference between the fine and 
coarse Greeks' estimates. This explains the experimental j3 w 1/2. 
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Figure 6: Vibrato Monte Carlo, digital call : V(Y/)(7) 



Value var with level Delta var with level 




3 4 5 6 7 8 3 4 5 6 7 8 



Vega var with level 




3 4 5 6 7 8 



4 European lookback call 

The lookback call's value depends on the values that the asset takes before expiry. Its payoff 

is P(T) = (S T - min (S,)). 

te[Q,T] 

As explained in [2], the natural discretisation P = (5jv — minS„) is not satisfactory. To 

n 

regain good convergence rates, we approximate the behaviour within each fine timestep 
[?„,?„ + i] of width hf as a simple Brownian motion with constant drift a n and volatility b{ 
conditional on the simulated values S n and S n+1 . As shown in [5] we can then simulate the 
local minimum 




with U n a uniform random variable on [0, 1]. We define the fine level's payoff this way 
choosing b{ r = b(Si,t n ) and considering the minimum over all timesteps to get the global 
minimum of the path. 

At the coarse level we still consider a simple Brownian motion on each timestep of 
width h c = 2hf. To get high strong convergence rates, we reuse the fine increments by 
defining a midpoint value for each step 

■W = \ {S c n +S c n+1 -K(AW n+l -AW n+l/2 )) (27) 
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Where (AW n +\ — AW n+ in) is the difference of the corresponding fine Brownian increments 
on [t n+ \/2,tn+i] and [t n ,t n +i/2]- Conditional on this value, we then define the minimum over 
the whole step as the minimum of the minimum over each half step, that is 



Sn,min ~ mm 



\ \S c n +S c n+l/2 - J(si +y2 -Stf-(btfh c logU lin ) , 



2 I 



J n+1 °m+1/2 



(bc) 2 h c logU Z ,n 



(28) 



where U\ in and U2 iH are the values we sampled to compute the minima of the corresponding 
timesteps at the fine level. Once again we use the tower property to check that condition 
Q is verified and that this coarse-level estimator is adequate. 

4.1 path wise sensitivities 

Using the treatment described above, we can then apply straighforward pathwise sensitivi- 
ties to compute the multilevel estimator. This gives the following results: 

Figure 7: Pathwise sensitivities, lookback call : V(Y/)(7) 



Value var with level 



Delta var with level 





Vega var with level 




Giles has proved that for the value's estimator, j8 = 2 — T] for all T] > 0. In the Black & 
Scholes model, we can prove that V = So 8. We therefore expected /3 2 for 8 too. The 
strong convergence speed of v's estimator cannot be derived that easily and will be analysed 
in our future research. 
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Table 7: Pathwise sensitivities, lookback call : estimated complexity 



Estimator 


P 


MLMC Complexity 


Value 


« 1.9 


<9(£- 2 ) 


Delta 


« 1.9 


0(£- 2 ) 


Vega 


« 1.3 


0(£- 2 ) 



4.2 Conditional Expectations, path splitting or Vibrato Monte 
Carlo 

Unlike the regular call option, the payoff of the lookback call is perfectly smooth and so 
therefore there is no benefit from using conditional expectations and associated methods. 



5 European barrier call 



Barrier options are contracts which are activated or deactivated when the underlying asset 
S reaches a certain barrier value B. We consider here the down-and-out call for which the 
payoff can be written as 

te[o,T] y 

Both the naive estimators and the approach used with the lookback call are unsatisfac- 
tory here: the discontinuity induced by the barrier results in a higher variance than before. 
Therefore we use the approach developed in [2] where we compute the probability p n that 
the minimum of the interpolant crosses the barrier within each timestep. This gives the 
conditional expectation of the payoff conditional on the Brownian increments of the fine 
path: 

N f -i 

pf = (S f Nf -K)+ J] (l-#9 00) 

n=0 

with 



p{ = ex P 



-2(S>-B)+(S' n+ l -B) 



At the coarse level we define the payoff similarly: we first simulate a midpoint value S° n+1 , 2 
as before and then define p c n the probability of not hitting B in [t n , t n+ \] , that is the probability 
of not hitting B in [t n ,t n+l / 2 ] and [t n+ y 2 ,t n +i]. Thus 

P c = (S c Nc -K) + if (1 -f n ) = (S c Nc -K) + if ((l-/5„,i)(l -p n ,2)) (31) 

n=0 n=0 

with 

-2(^-B)+(^ +1/2 -B)+\ 



p n>1 = exp 

Pn t2 = exp 



(K) 2 hf J 
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5.1 Path wise sensitivities 

The multilevel estimators % = (P^) ^ — (P c ) ^ ^ are Lipschitz with respect to all (S£) re =i..jv 
and (5^) n= i...Ar c , so we can use pathwise sensitivities to compute the Greeks. Our numerical 
simulations give 

Figure 8: Pathwise sensitivities, barrier call : ¥(?/)(/) 



Value var with level Delta var with level 




Vega var with level 




3 4 5 6 7 8 



Table 8: Pathwise sensitivities, barrier call : estimated complexity 



Estimator 


J8 


MLMC Complexity 


Value 


« 1.6 


0(£- 2 ) 


Delta 


^0.6 


0(£- 2A ) 


Vega 


^0.6 


0(£- 2A ) 
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Giles proved /3 = | — T] (T] > 0) for the value's estimator. We are currently working on 
a numerical analysis supporting the observed convergence rates for the Greeks. 



5.2 Conditional Expectations 



The low convergence rates observed in the previous section come from from both the dis- 
continuity at the barrier and from the lack of smoothness of the call around K. To address 
the latter, we can use the techniques described in section [T] Since path splitting and Vi- 
brato Monte Carlo offer rates that are at best equal to those of conditional expectations, we 
implement conditional expectations to see the maximum benefits we can get. 

Computing conditional expectations is slightly trickier than in section [2] We must in- 
deed take into account the probability that the path will hit the barrier B during the final 



timestep. Reusing the notations of part 2.2 and defining 



' a = 2B + (-\ + rh)S N -i(W) 
L = max(K,B) 

2r{B-S N -i) 
a 2 Vi 



D(o,S N -i) = exp 



(32) 



at the fine level we get 



>2k 



2[5 2 



(33) 



As before we then adapt this formula to the coarse level to compute 

E(P(S Nc )\W c ) = E(M(P(S Nc )\W c ,AW Nf -i)\W c ) 



(L-a 



rv\2 



(34) 



Doing so actually leads to long impractical formulae, especially when computing the Greeks. 
The idea of the conditional expectation technique is to smoothen the payoff. We can quickly 
estimate the method's maximum benefits by replacing the true payoff by a smooth Lipschitz 
approximation: this introduces a bias but also eliminates all the problems due to the lack of 
regularity around the strike. 



For example we can replace the payoff P = (St — K) 1 
approximation 



min (S t ) > 



g by the smooth 



'2% 



exp 



(K-S T ) : 
2p 2 



+ (S T -K)<S> 



S t -K 



1 min (5,) > K 
te[o.T] 



(35) 



where j8 = <jVJi*St for some arbitrary h* that controls the width of the smoothing. For 
example we take h* = 1/64 and we obtain the following results: 

We see in figure [9] and table [9] that the maximum benefits of these techniques are only 
marginal. The barrier appears to be responsible for most of the variance of the multilevel 
estimators. 
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Figure 9: Pathwise sensitivities and payoff smoothing, 
barrier call: ¥(?/)(/) 




Table 9: Pathwise sensitivities and payoff smoothing, 
barrier call : estimated complexity 



Estimator 


J8 


MLMC Complexity 


Value 


« 1.7 


<9(£- 2 ) 


Delta 


^0.7 


0(£" 2 - 3 ) 


Vega 


^0.7 


0(£- 2 - 3 ) 



5.3 Non-constant timestepping 



As illustrated in figure [10| the level at which V(F/) reaches its asymptotic convergence 
speed depends on the value of B. When B is far from So, the regime appears quickly (figure 



10a I, when B gets closer to So, it takes longer (figure 10b I. Practically this can be a problem 
when B « So as the simulations may not reach the very fine levels at which the complexity 
analysis based on the asymptotic value of /3 is relevant. 



In the case represented in figure 10b| the variance first increases before eventually con- 



verging towards 0. This illustrates the fact that in some cases it may be interesting to start 
the multilevel algorithm at a level / > where problems related specifically to coarseness 
do not appear and where the variance is low enough for our application (it must be at least 
lower than than the variance of the equivalent monolevel estimator). 

The variance's bad behaviour is related to the distribution of paths leaking out of the 
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Figure 10: Pathwise sensitivities (Vega), barrier call : V(F/)(/) 
(a)B = 85,5 = 100 (b) B = 95, So = 100 




barrier over time: we plot in figure 1 1 the density of first barrier crossings on the time 
interval [0,T]. We can show analytically that the width of the observed "peak" of the 



Figure 11: Density of first barrier crossings on [0,1] 



(a)B = 85,S = 100 



(b)B = 95,S = 100 



Distribution of 1 st crossing time, B- 85 



Distribution of 1 st crossing time, B= 95 




0.1 0,2 0.3 0.4 0.5 0.6 0,7 0.B 0,3 1 



0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,0 0,3 1 



crossing density function is 

t=0 (h*(y»)y (36) 

This means that as B tends to So, almost all paths going "down and out" do so in a very short 
time interval [0, t]. On timesteps outside of this interval, most paths are far away 

from the barrier and both ^ and ^ are close to and hardly contribute to the variance 
of the multilevel estimator. 

Morally the problem at the low levels is that the timesteps are much too large compared 
to the characteristic time v. the interval that is responsible for most of V(F/) is covered by 
only one step as long as hi > z. 

We hope to address this issue with adapted timestepping. Instead of taking constant 
timesteps of width h = T /N we use power timesteps to refine the discretisation in the time 
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interval [0, x]. We write t = u? and then split [0, r' 1 /^] into equal steps of u. We want more 
steps in [0, t]. Taking half of all timesteps in [0, t] means that we must choose y such that 
u = l/2 corresponds to z, that is 



iog(s /fl; 



log 2 



[logCT-log(log(5 /fl))] 



(37) 



Figure 
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shows for B = 95 the evolution of V(F,) for 8 (fig. [l2aj) and for v (fig. |l2b| 
with constant timesteps (red) and with power timesteps (blue). We see that with power 
timesteps, the variance is lower at fine levels and reaches its asymptotic convergence speed 
faster than before. Nevertheless we note that these benefits may be practically cancelled 
by the higher variance of the method at the coarsest levels, especially if we need rough 
estimates and stay at very low levels. 



Figure 12: V(7/), B = 95, power (blue) and constant (red) steps 



(a) 8 



Delta var with level 




(b)v 



Vega var with level 




Conclusion and future work 

In this paper we have shown for a range of cases how multilevel techniques could be used 
to reduce the computational complexity of Monte Carlo Greeks. 

Smoothing a Lipschitz payoff with conditional expectations reduces the complexity 
down to 0(e~ 2 ). From this technique we derive the Path splitting and Vibrato methods: 
they offer the same efficiency and avoid intricate integral computations. Payoff smoothing 
and Vibrato also enable us to extend the computation of Greeks to discontinuous payoffs 
where the pathwise sensitivity approach is not applicable. Numerical evidence shows that 
with well-constructed estimators these techniques provide computational savings even with 
exotic payoffs. 

So far we have mostly relied on numerical estimates of /3 to estimate the complexity of 
the algorithms. Our current analysis is somewhat crude ; this is why our current research 
now focuses on getting a rigorous numerical analysis of the algorithms' complexity. 
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